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Abstract 



We consider the infrared sensitivity of the inclusive heavy quark decay width in 
perturbation theory. It is shown by explicit calculations to the second loop order 
(when the non-abelian nature of the QCD interactions first become apparent) that 
all infrared sensitive terms to logarithmic and linear accuracy cancel when the decay 
width is expressed in terms of the short distance mass. This in turn implies the absense 
of any non perturbative I/uiq power corrections to the decay width. The cancellation 
is shown algebraically at the level of the Feynman integrands without any use of an 
explicit infrared cutoff. This result is in accord with the implications of the KLN 
theorem. 



1 



I. INTRODUCTION 



The theory of heavy quark decays has attracted much interest lately in connection with 
the decays of B mesons containing a single b-quark. Of particular phenomenological interest 
are the inclusive decays which can provide information about the CKM parameters. Inclusive 
decays of hadrons containing a single heavy quark (mass rrib) have been extensively studied 
1 1 J] within the operator product expansion (OPE) which gives an expansion of the decay 
width in terms of a small parameter A/m& : 

T ind = rf n e c r [(l + a 1 A/m 6 + a 2 (A/m & ) 2 + ...). (1) 

In the above, T^ r J is the decay width of the heavy quark including all perturbative correc- 
tions (the parton model result). Since there are no local gauge invariant operators of the 
appropriate dimension, an important conclusion Q is that ai = 0, and the higher order 
power corrections can be classified in terms of the matrix elements of various operators. 

A great advantage of the OPE is the generality of its applicability to non-perturbative 
contributions as well. On the other hand, the conditions for the validity of the OPE itself 
can be enunciated by explicit perturbative calculations. As far as the power corrections are 
concerned, the method of renormalons (for a review and further references see ||) may be 
used to study these within perturbation theory. The basic idea of the renormalon method is 
that one may obtain information about power corrections by looking at a class of diagrams 
which give factorial divergence in large orders of perturbation theory. The power corrections 
are seen to arise from regions of low momenta and thus we are led to investigate the infrared 
sensitivity of inclusive decays. To this end A, which is a non-perturbative parameter is 
replaced by an infrared cutoff defined within perturbation theory. We will generically denote 
this infrared cutoff by A and in the case of QED, for example, it could be a (fictitious) photon 
mass. Then the inclusive rate ([!]) may be represented as : 

Find = TZ c \(l + b l \/m b + b 2 (\/m b ) 2 ln\ + ...). (2) 

Note that the Landau conditions for the singularities of Feynman diagrams tell us that the 
infrared sensitive contributions arise only as terms which are non- analytic in A 2 . Explicit 
calculations || to leading order in a have brought to light a subtlety that the coefficient 
bi is zero in perturbation theory only if the decay rate on the right hand side of Eq. @ 
is expressed not in terms of the pole mass m& but rather in terms of a short distance or a 
MS mass. Thus perturbation theory is in full accord with the OPE, only if in ([I]) the short 
distance mass is used. In this way one does not have to worry about an otherwise large 
(possibly non- universal) correction of order 1/mt,, but only at the expense of introducing a 
renormalization scale (/i) dependence in the total rate through the use of the running mass. 
We will discuss below that this result arises from the cancellation of the infrared sensitive 
contributions proportional to A/m& (and of course to InX) arising from the long distance 
pieces of the mass and wave function renormalizations of the initial heavy quark and from 
the gluon radiation. 

The result mentioned at the end of the previous paragraph has a simple physical inter- 
pretation. The decay rate may have two possible sources @] of infrared sensitivity : (a) the 
heavy quark mass and (b) the energy radiated during the decay. To expose the infrared 
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sensitivity of the heavy quark mass say within QED, consider a QQ pair separated by a 
distance r and interacting by the exchange of photons of mass A. The total energy is 

E = 2m b - — exp(-Xr) ~ 2m 6 (l + < ^1—) - a d /r (3) 

where we have kept only the term linear in A. It is clear from this that we may identify 

m[ R = a el X/2, (4) 

with the infrared sensitive piece of the heavy quark mass and for small r it is appropriate 
to define the short distance mass rh by: 

m = m b (l + (5) 

ZlTLb 

The same result was obtained in [[3] by the usual renormalon method . Let us next consider 
the radiated energy. If the decay time is r then classically, a light signal cannot reach a 
distance beyond 

Rcausal = C ■ T (6) 

As a result, the field energy stored at distances larger than R caU sai cannot be transformed 
into the energy of the decay products and cannot affect the decay rate and hence must be 
radiated away. Let us denote this long distance piece of the energy radiated during the 
decay by E^ d . That this long distance part of the radiated energy is proportional linearly 
to the infrared cutoff may be seen from the classical expression for it: 

EZ = f dul(uj). (7) 

J small u> 

I(uj) is the intensity of radiation which approaches a constant || as uj — > 0. The cancellation 
of the linear term in Eq. (0) is the statement that to this order, 

m{ R = E™ (8) 

and the decay rate which is sensitive to the energy release at short distances is not affected. 

Certain questions immediately come to mind concerning the generality of these results. Is 
there a general principle behind this cancellation ?, and is it true to all orders in perturbation 
theory, even for the non-abelian case? In fact there are infrared safe observables in QCD that 
do receive power corrections that are proportional to 1/M where M is a large mass scale. 
Examples are provided by the event shape variables in e + e~ annihilation for which there 
does not exist an operator product expansion. It has been argued in ref. [0] that inclusive 
enough observables do not receive power corrections of the type 1/M whereas those that 
assume some precision measurement (and hence are more exclusive) like the event shapes, in 
fact do. Below we will apply this general principle to the inclusive decay of the heavy quark. 
We will argue that the KLN || theorem is behind the cancellation of the terms proportional 
to l/mjj. The purpose of this paper is to show by an explicit calculation in the full QCD 
theory the cancellation of the infrared sensitive pieces that are linear in the infrared cutoff 
to order a 2 in the decay rate. 
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In a recent publication || it was argued that when the transition probability is summed 
over both initial and final degenerate states as is envisaged in the KLN theorem, not only 
the leading infrared sensitive piece proportional to InX but the one proportional to A is 
cancelled as well, i.e, 

\Si->f\ 2 ~ • InX + • A + terms independent of A + 0(X 2 lnX). (9) 

i,f 

The transition probability summed over both initial and final states cannot be related to 
a physical inclusive cross section. However we will now argue that for heavy quark decay, 
the initial state is trivial and hence, the above mentioned property of the KLN theorem 
guarantees the cancellation not only of the In A terms but also those proportional to A itself. 
To substantiate this we note that the KLN theorem presumes that the total energy can be 
fixed to any accuracy. This is in the spirit of the uncertainity principle, which allows us to 
measure the total energy to any accuracy provided the measurement time is long enough. 
In particular, the uncertainity AE in the total energy can be made smaller than an infrared 
cutoff: 

AE < A (10) 

Now consider a decaying charged (or colored) particle in its rest frame. It is obvious that 
this state is not degenerate with any other since there is a gap between the mass of the 
charged particle and the energies in the continuum for any A ^ 0. This is most clearly seen 
for QED with a non-zero photon mass, but is true for any infrared cutoff on the energy, A. 
Thus the summation over the initial states in the KLN theorem is redundant for this case. 
The final states however include the full degeneracy but the summation over the final states 
is equivalent to calculating the decay width, and hence we conclude that, 

r,nd = TZl(l + 0(\ 2 ln\)). (11) 

It should be emphasized that this argument guarantees the above equality to all orders in 
perturbation theory. However we should note that it is very important for the KLN theorem 
to be valid that the renormalization procedure does not introduce any infrared sensitivity. 
Thus since the pole mass is infrared sensitive, on shell mass renormalization must be avoided 
in favour of the MS scheme . In concluding this we would like to point out that the same 
argument can be applied not just to the semileptonic decay but to other inclusive decays 
like for example the radiative one, B — > X s 7. In fact, the explicit perturbative calculation 
though performed for the semileptonic decay, can equally well be applicable for the radiative 
decay b — > s^y with minor changes for the electroweak vertices and phase space factors. 

As mentioned previously, we will explicitly verify the above argument in this paper for 
the inclusive semileptonic decay of a heavy quark in QCD to the second loop order when the 
non-abelian nature of the interactions are first apparent. That is we show the cancellations of 
all terms linear in the infrared cutoff in the decay width when the latter is expressed in terms 
of the short distance mass. The possible sources of infrared sensitive terms are the mass and 
wave function renormalization diagrams of the initial heavy quark and the bremsstrahlung 
diagrams. We will discuss the cancellation of the logarithmic and the linearly infrared 
sensitive terms in the diagrams for the wave function renormalization and of the linearly 
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IR sensitive terms in those for the mass renormalization. These cancellations proceed in 
a different manner for the two cases because whereas the wave function renormalization 
is multiplicative, the on-mass shell renormalization is additive. In section [Tl] we introduce 
our notations and conventions and discuss our strategy. In particular we show how to the 
linear accuracy desired, we may replace the diagrams for the radiation from the final state 
massless quark by effective local vertices. Next we discuss the wave function renormalization 
in section |TT| and the mass renormalization in section [TVJ Certain additional technical details 
are relegated to the appendix. In section [V| we summarize our results and discuss future 
prospects. 



II. PRELIMINARIES 

We consider the totally inclusive semileptonic decay of a heavy quark and examine the 
possible term linearly proportional to an infrared cutoff. Such a term can arise from the soft 
gluon contributions to the diagrams for the mass and wave function renormalizations of the 
heavy quark and from the bremsstrahlung diagrams. 

The decay rate is proportional to the imaginary part of the forward amplitude shown in 
Fig. 0. The hard amplitude involving the final state quark is denoted by the shaded blob 
which also includes the lepton loop. The soft gluon interactions responsible for the infrared 
sensitivity dress the hard amplitude, connecting to it and to the heavy quark. We will use 
the fact that for the soft gluon of momenta k <C m^, we may perform the sum over cuts 
implicit in taking the imaginary part, for just the hard part of the amplitude leaving the soft 
gluon lines uncut. This simplification arises due to the fact that for the process of a heavy 
quark decay, we may assign a positive energy flow direction to each soft gluon consistently. 
Then we may write for the corresponding gluon propagator: 

2n9(k )5(k 2 ) + - ,V (12) 



k 2 + ie k 2 + ikoe 

The second term can never cause a pinch singular point since the contour can be deformed 
into the lower half fco plane. This piece therefore (from the Landau equations) does not give 
any contributions that are non-analytic in the infrared cutoff. Thus as far as the infrared 
sensitive terms are concerned, even upto linear accuracy it does not matter if we cut the gluon 
propagator or not, both possibilities giving the same contribution. Keeping this in mind, 
we will sum over the cuts of the hard amplitude only and leave the soft gluon propagators 
uncut in the diagrams for the forward scattering amplitude. The next simplification of the 
forward scattering amplitude we will use for our infrared analysis is that the (heavy quark) 
propagators on the right hand side of the cut can also be taken to have the +ie prescription, 
the same as those on the left hand side of the cut. The difference between the two choices 
is readily seen to be pure imaginary and since the cut hard amplitude is already imaginary, 
we may ignore this difference for our purposes. 

Let us consider the diagrams with the soft gluon radiation interacting with the hard 
part. In the next subsection we will show that to the accuracy we want, after summing 
over the cuts of the hard part, and integrating over its phase space and that of the leptons, 
this interaction may be replaced by effective local vertices for the interaction of one and 
two gluons. These effective vertices may be obtained from the corresponding one without 
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gluon radiation by gauge invariance. Thus we can forget about the details of the hard 
part and use the effective vertices to discuss the infrared sensitive contributions for such 
diagrams. Diagrams in which all the gluon interactions are contained in the hard part 
are irrelevant to our analysis as far as the infrared sensitive terms are concerned which 
are linearly proportional to the infrared cutoff. Since the wave function renormalization is 
multiplicative, all its contributions factorize from the hard amplitude. In addition, as we 
discuss below the corresponding bremsstrahlung diagrams that cancel the logarithmic and 
linear dependence on the IR cutoff of the wave function renormalization diagrams may be 
obtained by certain Ward-like identities. These cancellations can be shown independent of 
the hard amplitude whose precise structure is not needed here, but the latter is essential for 
showing the cancellations among the mass renomalization and the corresponding radiation 
type diagrams. The use of the heavy quark rest frame provides very important simplifications 
of this analysis. 

In this paper all the infrared cancellations are shown strictly algebraically at the level 
of the corresponding Feynman integrands. Power counting is used to isolate and to show 
the cancellations of the possibly infrared sensitive contributions upto linear accuracy. After 
the cancellations, from power counting, the integrand is seen to contribute something which 
would at least be proportional quadratically to the infrared cutoff. In this way we avoid 
using any particular IR cutoff. This is crucial since for the non-abelian theory we do not 
know of any gauge invariant infrared cutoff for explicit calculations at the two loop level. 

Finally we note that throughout this paper we work in the Feynman gauge. The can- 
cellation of the infrared sensitive terms is a gauge invariant statement, however, the type of 
diagrams contributing to this cancellation depends on the choice of gauge. 

A. Effective vertices for the Interaction of radiation with the Hard Amplitude 

In this section we will derive the effective vertices for the interaction of the gluon radiation 
with the final state massless quark in the hard amplitude.. These effective local vertices are 
constructed to reproduce the soft momentum (k <C m^) contributions upto linear in l/m b 
accuracy. Consider the absorbtion of a single gluon of momentum k from the final state 
quark. We will show below that summing over the cuts for the final state fermion line cancels 
the term proportional to 1/k (i.e., the would be infrared divergent term proportional to ln\ 
) so that the leading left over piece is of order k°. When the other end of the gluon line is 
attached to the heavy quark line then we would be interested in the term proportional to 
A and from simple power counting of the corresponding diagram we see that we need only 
keep this leading term proportional to k° from the sum over cuts of the final state quark. 
In this way we will conclude that after the phase space integration over the leptons and the 
final state quark at fixed k this effective interaction may be represented by a local vertex 
of the form shown in Fig. |3|b). We show below that this effective local vertex is what we 
would obtain by the substitution p — > (p — gA) in the leading effective vertex of Fig. |3]a). A 
similar effective vertex is derived for the absorbtion of two gluons from the final state quark. 
Since as mentioned earlier, the heavy quark wave function renormalization contribution to 
the infrared sensitivity is multiplicative in nature, these contributions are not sensitive to 
the details of the interaction of the final state quark and the effective vertices derived in this 
section are not needed there. 
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In the parton model, the total width for the semileptonic decay is to lowest order: 



1927T 3 



■mi 



(13) 



We may represent this by an effective local vertex, obtained by preforming the phase space 
integration over the leptons and and the massless quark by the heavy dot in Fig. |3|a). Such 
an effective vertex may be written as: 

Cufu (14) 

where, u is the heavy quark spinor. In the above, 

r 2 IT/ |2 

c - ^L. (15) 

The vertex of Fig. [3|a) will be referred to as the leading effective vertex. Next consider the 
class of diagrams of Fig. 0a). The gluon may undergo self interactions before connecting 
to the heavy quark line. Thus, for example, the diagrams of Fig. [|b) are included in the 
set. Next we sum over the two cuts on the massless quark line. In one just the quark line is 
cut and in the other, both the quark and gluon lines are cut. However, as discussed above 
and as it was shown in explicit examples in ref. 0, as far as the infrared sensitive terms 
(those that are non-analytic in the infrared cutoff) are concerned, both real and virtual soft 
quanta may be regarded the same i.e, we may write for the gluon line contribution either as 
2n9(k )5(k 2 ) or if (k 2 + ie). This is really a consequence of the Landau equations as discussed 
earlier. Keeping this in mind we can write for the sum over cuts in Fig. |2|a) as: 



-ry 11 a g T a (j/ 1 -^)r27r 



5{p\) , 5{{ Pl - kf 



(pi — k) 2 + ie pi — ie 



(16) 



where T generically denote the weak interaction vertices, and p\ = (p — q) with q the 
momentum transfer to the leptons. Note that the lepton loop is always cut and we do 
not explicitly write down the corresponding expression. We also have not explicitly written 
down the contributions from the blob in Fig. |2]a) and those from the external quark lines 
since their exact form is inessential to the argument. Now we can use the following identity 
to simplify the term in squared brackets above: 

5{p\) + 5(( Pl -k) 2 ) 



(pi — k) 2 + ie p 2 — ie 
(—1) J dx\ J dx 2 S(l — Xi — x 2 )S^ [x\p\ + X2(pi — k) 2 ^j . (17) 

In the limit k — > the right hand side of Eq. (|i~7l) becomes: 

(18) 

which is a well defined quantity. This shows the cancellation of the leading infrared divergent 
term, and the next term in the sum over cuts, which is 0(k°) may be written as: 
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-rVi7 Q ^Vir2vr(-l)5( 1 )(p?) 



(19) 



We wish to compare this with the corresponding expression for the absorbtion of a zero 
momentum gluon from the uncut final state massless quark of momentum pi which is: 



Using the fact that: 



^WrW' (20) 



Mrn = 2^ 1 \pj), (21) 

(p( + ie) z 



and comparing with Eq. (|19D we see that as far as the terms of order k° are concerned we 
may represent the diagram of Fig. 0a) for fixed k as an effective local interaction which 
after the appropriate phase space integrations over the leptons and the massless quark for 
a fixed k may be obtained from the effective vertex Eq. fll4|) by means of the replacement 
P/j, —* (Pfj, — gT a A a ^). We have of course shown this only for the one gluon effective vertex 
shown in Fig. |]b). Next we will show that the same substitution also reproduces the two 
gluon emission or absorbtion process. 

Consider next the absorbtion of two gluons by the final state quark as shown in Fig. 0c) 
and once again the blob represents all possible interactions and connections of these gluons 
with each other and with the heavy quark. The sum over cuts of the massless quark line for 
this configuration for fixed fci, k 2 is: 



r^x^v (#! - fa) gTV {ii -fa- fa) T27T 
i i 



(Pi - k 2 ) 2 + ie (pi —ki — k 2 ) 2 + ie 

S ((pi - k 2 f) -^—- ( 1 . + 5 ((p! - h - k 2 f) 7ry 2 : 

v ' p{ — ie (pi — ki — k 2 ) A + ie v ' p{ — ie [pi — k 2 y — %' 

Next we use a generalization of the identity Eq. (0), 

11 11 

^l)? ? \ 2 I • 7 ; 1 \ 2 I + 6 ((Pi - k l) 2 ) — -7 j 1 \ 2 I + 

(Pi - k 2 ) z + ie (pi — Ki — k 2 ) + ie v y Pi - ie (pi — ki — k 2 y + ie 

8 {(Pi - h - k 2 ) 2 ) -= — = (-1) 2 / .. / dxidx 2 dx 3 5 (1 - x x - x 2 - ar 3 ) x 

v y pf — 2e (pi — K2) — «e Jo jo 

^ (xxp 2 + x 2 (p x - £; 2 ) 2 + x 3 (pi - h - k 2 ) 2 ) (23) 



(22) 



The right hand side of (p23|) tends to a well defined limit as ki —> 0, 

(-1) 2 ^ (2) (P?) (24) 

which again represents the cancellation of the leading logarithmically divergent contribution. 
The term proportional to k° in the sum over cuts is then given by: 

Y^ 1 gT a ^ 1 gT b 1 ^ 1 T{-l) 2 ^ 2 \p 2 l ). (25) 



S 



Again consider the corresponding term for the absorbtion of two zero momentum gluons 
from the massless quark line: 

- i^h9T a l a ^ 9 TV^ {pl l ief (26) 

Since 

2Um MTw = - n6(2) ^ (27) 

comparing with ([25|) we see that the 0(k°) contribution from the sum over cuts may again 
be written as an effective local vertex after the relevant phase space integrations. This local 
vertex can be obtained from Eq. (|14D by means of the replacement p M — > (p^ — gT a A a ^) and 
expanding to 0(g 2 ). Note that there is a similar contribution as in Eq.(|25|) but with the 
interchange (a, a) <-> ((3, b), and this too is reproduced by the above mentioned replacement. 
The two gluon effective vertex is shown in Fig. |3]c). 

We should remark here that when we go beyond the two gluon emission case, in general 
multi-gluon emission or absorbtion from the final state line cannot be written as an effective 
local vertex. Fortunately this is not a problem for us. 



III. CANCELLATION OF INFRARED SENSITIVE TERMS UPTO LINEAR 
ACCURACY FROM THE WAVE FUNCTION RENORMALIZATION 

DIAGRAMS 

In this section we will show that the infrared sensitive terms upto linear accuracy are 
cancelled between the wave function renormalization contributions from diagrams with self 
energy insertions and parts of other diagrams related to the bremsstrahlung. The cancelling 
sets can be arranged into groups which as we will see below upto linear accuracy are obtained 
through certain Ward-like identities. For the cancellations considered in this section which 
are related to the wave function renormalization, the corresponding contributions are such 
that the hard part of the diagrams factorize. The hard part thus is not consequential to the 
argument and just comes along for the ride. 

Throughout the paper, we will denote by — the sum of all 2 point one particle 
irreducible graphs. Then the full fermion propagator is: 

f -m b - E(p) 

Near the mass shell we have 

S(p) = 6m (m b , Am) - (Z' 1 - 1) {p 1 - m b ) + O ({p 1 - m h f) (29) 
where Am denotes the mass counterterm: 

5m (m b , Am) ^ = 0. (30) 
Thus m b denotes the pole or the physical mass: 
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lim S' F (p) = Z 

■Hm b f — nib 



(31) 



From the LSZ reduction formula, each external line in the S-matrix element gets a factor 
Z 1 ! 2 and we are interested in its perturbative expansion. Thus we expand, 

Z 1 ' 2 = l- 1 -(z^-l) + 3 -(z-i-l) 2 + (32) 

and Z^ 1 — 1 may be computed using: 

Z~ x -\ = - dVv ^ mb (33) 
wy v u 

In the rest frame of the heavy quark, only v — contributes. 

We will now clarify our procedure first with the 0(a) case and next the two loop example 
will be discussed. 



A. Cancellation to One-Loop Order 

Consider the graphs of Fig. f|. The wave function renormalization parts of both graph(a) 
and graph(b) are the hard part times a factor of Z 1 ! 2 . Thus we expand it using (|32|) keeping 
only the first two terms and computing the perturbative contribution to order a of Z~ x — 1 
from (|33|). Next using the identity 

d i i . i 

Qpv ^ _ ^ _ mfe ^ _ ]j- _ mb — — TTlb 

we obtain the Ward-like identity shown in Fig. |, where the dotted line just denotes the 
insertion of a zero momentum vertex r y u . We would now like to show that upto logarithmic 
and linear accuracy in the infrared, the right hand side of Fig. ||] may be identified with the 
wave function renormalization part of Fig. |4]c) times u^ v u, modulo the hard part, which 
after the relevant phase space integrations just gives a factor of Cm\. It is easiest to go now 
to the rest frame of the heavy quark, where 

7°m = u. (35) 

From this and Eqns. (0), (0) we immediately see 

4c) \ wavefuncUo n = (Z- 1 - \)?Cm\uu = - 10 + BO] ■ (36) 

Here, and henceforth, unless otherwise specified, the equality sign is for terms that are loga- 
rithmically and linearly infrared sensitive. This is the desired result, which does not involve 
using an explicit infrared cutoff. Let us now see how the wave function renormalization part 
of Fig. [|c) can be identified with (Z^ 1 — l^Cmjj, which is effectively the right hand side 
of Fig. |5| multiplied by the hard part Cmj. Consider the expressions for these diagrams in 
the rest frame of the heavy quark. Both are of the form: 
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{ig) 2 C F Cml { -f^-N % —— % — (37) 

v yj b J (2vr) 4 -2pk + k 2 -2pk + k 2 k 2 K ' 

where N denotes the numerators N$ or N± c for the two diagrams respectively. Thus we only 
have to compare the numerators to linear accuracy. Upto linear terms in the infrared power 
counting 

N 5 = W7 a (#f + m)7°($f + m)^ a u - u^ a ^(^> + m)^ a u (38) 
—u^ a {jj> + m)7°^7 Q ,M 

There is a similar expression for the numerator N 4c except for the factor of 70. However, 
using Eq. ( |35l) it is easy to see that to this accuracy the two numerators are equal 

N 5 = N 4c (39) 

This demonstrates the cancellation of not only the infrared divergent part but also the next 1 
inearly infrared sensitive piece in the wave function renormalization. After the cancellation, 
the integral will be at most quadratically divergent by IR power counting. The cancellation 
of the wave function parts at two loops follows the same general method except that now 
we have more diagrams. 



B. Cancellation to Two Loop Order 

The 2-loop one particle irreducible diagrams can be divided into 5 groups according to 
the 5 pieces of the order g A self-energy corrections (Fig. |7|). These groups are essentially 
the different color groups, with the exception that the first and third group have the same 
color structure. In this section we show the cancellation of infrared sensitive terms to linear 
accuracy in the second loop order for the wave function renormalization diagrams and their 
radiation counterparts. Power counting in the IR region is used to isolate logarithmically 
and linearly infrared sensitive pieces of the integrals. In the following, these infrared sensitive 
terms identified by power counting will be simply referred as "logarithmic and linear pieces" 
of the corresponding diagrams. As it is discussed below, the cancellation of logarithmic and 
linear terms takes place within each color group separately. The cancellation is algebraic at 
the level of integrands and no explicit infrared cutoff is used. 

Let us take the second group (Fig. ||) as an example to work with. We would like to 
show that both the logarithmic and the linear pieces in the wave function renormalization 
of diagrams |9]a) and ||b) cancel the corresponding logarithmic and linear wave function 
renormalization pieces of diagrams |9|c), |9]d) and []e), i.e. 

5 ( rv-\ 



[9a) + miavefuncUon = ~ Cm b " «« (40) 

= — @c) + ^p) + ^f)] wa vefunction 



As before, we may write 



u 



1 



u u 

i>=m b 



22 



9 4 .£22 wy u u uu 



u 

i^m b _ (41) 
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where the last equality holds in the rest frame of the heavy quark p = (nib, 0, 0, 0). 
We would like to show that to linear accuracy 



Cmiu 



5 _ 8T.22 



dp 



Since 



u = iCm\u 

jt=m b 

-M+W)+W)] W avefunction 







_ d(-iZ 22 ) 



dp 



u 



(42) 



jj=m b 



- ii° - (43) 

dpo i> — k 1 — m f — If — m f — Jp — m 

the derivation of £22 with respect to po results in the insertion of an extra vertex ry° and 
an extra propagator in every possible way. The resulting diagrams are similar to the wave 
function renormalization parts of |9]c), |9|d) and|9|e) with the Cm\ vertex replaced by — Cm^7° 
(see Fig. |l~3|b). Writing out explicitly 

CpCa 



5 _ <9£ 2 2 

Cm b u 

"0 



dp, 



u 



-Cm\ C 2 



t>=m b 



-) (ig) 4 x I J d 4 kd%N r 



-2pk + k 2 



(44) 



—1 —1 



-2pk + k 2 -2p (k + h) + {k + h f -2ph + k\ k 2 k 2 

i i i 



+ / d 4 kd%N 2 x 



—1 —1 



-2pk + k 2 -2p (k + h) + (k + A7) 2 -2p (k + h) + (k + A7) 2 -2ph + k\ k 2 k\ 

f 1 % % % 

+ / d 4 kd 4 k\N-i — ^ -77 ^ 

J -2pk + k 2 _2p(fc + jfc 1 ) + (jfc + fc 1 ) 2 -2ph + k 2 -2ph + k 2 



-1 —1 



k 2 k\ 

where Ni, N 2 , N 3 are the numerator factors 

Ni = u-f a {J - k 1 + m) 7 - k 1 + m) 7^ {p 1 - k 1 - fa + m) -? a (p 1 - fa + m) lp u (45) 
N 2 = u^ a (p 1 - Jjt + m) 7^ (#f - Jjt - fa + m) 7 (p 1 - Jjt - fa + m) 7 Q (p 1 - fa + m) jpu 
N 3 = u^ a - Jjt + m) j p (J - - fa + m) j a {$ - fa + m) 7 {p 1 - fa + m) j^u 

Comparing this with diagrams [|c) , |]d) and ||e) 

CpC 'a \ / . \ A f ,A , ,A , t> t i 



-9c) 



-Cml ( C| - 



{igf J d 4 kd 4 k x N Xc - 



2pk + k 2 
—i —i 



(46) 



-2pk + k 2 -2p (k + h) + (k + hf -2ph + k\ k 2 k\ 
-BO = -Cm\ (Cl - ^) (^) 4 / d 4 kd%N 2d 1 



x 



2pk + k 2 

i —i —% 



-2p (k + h) + (k + fa) 1 -2p (k + h) + (k + h) 2 -2ph + k\ k 2 k\ 
-|e) = ~Cm\ (C 2 F - (igf [ d 4 kd%N 3e - ' 



-2pk + k 2 
-i —i 



x 



-2p (k + fci) + {k + h) 2 -2ph + k\ -2pk x + k\ k 2 k\ 
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with the numerator factors 

iVic = ~ $ + m) - lfi + m) ^ (j> - # - fa + m) la [j> - fa + m) 7/9 u (47) 

N 2 d = u^ a {j> - Ijt + m) 7^ {j> - $ - fa + m) (J - # - fa + m) ^ a {j> - fa + m) 'jpu 
N 3e = u>y a {i> - Jjt + m) 7^ (|f - # - fa + m) 7 Q (jl - fa + m) - fa + m) 7/3 u 

we see that the proof of equation ( f4"2"D is tantamount to showing the equality of numerator 
factors to linear accuracy. IR power counting on the integral shows that we have to keep 
the terms in the numerator at most linear in the integration variables. To linear accuracy 
the first numerator pieces, N\ and N\ c are 

Ni = M7 a (J + m) 7° (J + m) 7^ (J + m) 7a {j> + m) ^pu (48) 

- M7 a ^7° (rf + m) 7^ {$ + m) 7 Q (^ + m) ^pu 

- «7 a (3/ + m) 7°^7 /3 + m) -y a (tf + m) ^pu 

- «7 a (jrf + m) 7° (jf + m) 7^ + ^) 7^ (^ + m) ^pu 

- «7 a (^ + m) 7 (^ + m) 7^ (^ + m) j a fajpu 

and 

^lc = M7 a (#f + m) (^ + m) 7^ (#f + m) 7^ (#f + m) ^pu (49) 

- tt 7 a # ($f + m) 7^ (^ + m) 7 Q (#f + m) ^pu 

- uj a + m) ^7^ (^ + m) 7 Q (#f + m) jpu 

- u^ a + m) + m) 7^ + ^1) 7q + m) ^pu 

- uj a (J + m) + m) 7^ (^ + m) ^atyi^pu 

Note that the first terms in equations ( fl8|) , ( ^9|) give the logarithmic infrared sensitive pieces, 
while the last 4 terms give the linear infrared sensitive parts. Comparing the first terms, it 
is easy to see the equality of the logarithmic pieces. 

It is also straightforward to see that Ni\ Unear = N lc \ Hnear . The only difference in the 
numerators is the presence of the 7 factor. Since the numerator is at most linear in k and 
hi, either left or right to 7 there will be only factors of + m) and 7 matrices, which 
could be easily evaluated using the mass-shell condition. Then using that 7°m = u (since in 
rest frame) we have effectively replaced the 7 by 1, which means the equality of the two 
numerator pieces. In the same way one can show that to linear accuracy N 2 = N 2 d and 
^3 — -^3e; thus the equality of numerators and hence equation ([42]) holds to linear accuracy. 
Therefore the logarithmic and linear infrared sensitive wave function renormalization pieces 
of the £ 2 2 group cancel 

[9a) + 16) + |C) + gd) + MwauefuncUon = (50) 

The cancellation of the logarithmic and linear wave function parts of group-4 can be 
shown in a similar way. The cancellation of logarithmic and linear wave function renormal- 
ization pieces in group-5 is discussed in the Appendix. 

The first and third group do not cancel separately, since S 2 i contains a mass renormal- 
ization piece through its inner part, Ex. This piece will cancel with the corresponding mass 
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counterterm of E 2 3, and the sum of the two groups again cancel to the accuracy desired. 
To show the cancellation of the logarithmic and linear infrared sensitive pieces in the wave 
function renormalization, we again use the Ward-like identities in Fig.s |i~3|a) and [13|c). 
At the two-loop level 1-particle reducible diagrams also contribute to the wave function 



renormalization. The relevant diagrams are shown in Fig. [14]. The contribution of the 
diagrams (shown in the Figure) can be simply calculated using the perturbative expansion 
of Z 1 ! 2 (Eq. |31| ) multiplying each external line in the S-matrix element. As Fig. [TJ]shows, the 
wave function renormalization pieces from the two-loop order 1-particle reducible diagrams 
cancel. 



IV. CANCELLATION OF THE INFRARED SENSITIVE TERMS ARISING 

FROM MASS RENORMALIZATION 

Let us now consider the cancellation of the leading infrared sensitive piece (the one 
that is linearly divergent in the infrared) arising from the diagrams involving the mass 
renormalization of the heavy quark. As for the wave function case we begin with a discussion 
of the one loop case to exemplify our method and then the more complicated two loop 
example is considered. 



A. Cancellation at One Loop 

As discussed in the introduction, the pole mass of the heavy quark contains a long 
distance piece which is proportional linearly to an infrared cutoff, the first contribution 
starting at order a. In fact, for our subsequent analysis it is convenient to separate out this 
piece and write for the pole mass to one loop order: 

m h = fh + big 2 \ (51) 

where, as before, fh is used to denote the short distance or the running mass, which itself is 
a power series expansion in g 2 : 

fh = m + dig 2 . (52) 

We should emphasize our notation at this point. Even though one may use a gluon mass 
as an infrared cutoff at the one loop level, this is not a gauge invariant procedure beyond 
it. Thus A refers to some gauge invariant cutoff. Throughout this paper, however we will 
not have any need to specify this infrared cutoff. Terms that are linearly divergent in the 
infrared are identified by power counting and cancellation of such terms is shown at the level 
of the corresponding Feynman integrands. Thus for example, 

b l9 2 X = g 2 C F [ ^t^—tt^- (53) 

Jlinear (/VTJ 4 K z + 16 K — It 

The various diagrams contributing to this order are shown in Fig. ^. Fig. |]a) is the 
bremsstrahlung contribution. To find the leading infrared sensitive contribution, we expand 
the leading effective hard vertex to first order in the offshellness, i.e., 
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C{i> - Iff = C{i> - # - m b + m b f = Cm\ + C5m A b (tf - # - m b ) + 



(54) 



This is an expansion of the leading effective vertex in {jj> — If, — m b ) / m b and it is understood 
that we keep such deviations from the mass shell in as much as they cancel the corresponding 
small denominators from the propagators. Then after this expansion, one of the fermion 
propagators in the integrand of the expression for Fig. ^ is cancelled (shown by a slash in 
the Figure) and we are left with an expression which resembles that for the self energy. In 
fact, 

6a) = C5(rh)%g 2 \uu (55) 

Next we have the contributions from the bremsstrahlung from the final state quark which 
as discussed earlier is replaced by the effective vertices for the gluon emission. The two 
diagrams with single gluon effective vertices of Figs. |5|b) and |^c) give: 

66) +|c) = 2 x C(-ihm A b )u (-i£i(p)) u = -2 x C5(rh)\g 2 \uu. (56) 

Finally, in the lowest order diagram, Fig. |6]d), we must express the leading effective vertex 
in terms of the short distance mass, 

Qd) = Cm\ = C(m + b ig 2 \f = Cm 5 uu + 5C(m)%g 2 \uu. (57) 

Adding together all the order g 2 terms we see the cancellation of the linearly infrared diver- 
gent pieces to this order when the decay rate is expressed in terms of the running mass rh. 
We next discuss the two loop example. 



B. Cancellation to Two Loops 

As in the wave function renormalization, we group the two-loop 1 PI diagrams according 
to the five self-energy pieces shown in Fig. [7]. These groups are now enlarged by new 
diagrams with the single and double gluon vertices (Fig. |^). Unlike in the wave function 
case, where each group of graphs cancelled separately, the cancellation of linear infrared 
sensitive pieces in the mass renormalization involves mixing between the different groups 
and mass shift terms from one loop graphs. This is essentially because whereas the wave 
function renormalization is multiplicative, the on-shell mass renormalization is additive. 

By infrared power counting we can isolate the linear infrared sensitive pieces of the 
self-energy diagrams as 

Si|^ mb =a ig 2 + b ig 2 X (58) 
s 2i|^ mfc = a 2i g* + b 2 ig A X i = 1 . . . 5 

(Note that as before, A denotes symbolically the linearly infrared sensitive pieces of the 
diagrams identified by power counting, not an actual infrared cutoff parameter.) Define the 
running mass as 

5 

fa = m + a x g 2 + ^ a 2i g 4 . (59) 

i=l 
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Then the physical mass in terms of the running mass can be written as 

5 

m b = m + b 1 g 2 X + b 2 ig 4 X. (60) 



8=1 



1. Group- 1 and Group-2 

a. Leading effective vertex and single gluon vertex graphs The cancellation between the 
linear infrared sensitive pieces in the mass renormalization of the first two group of graphs 
with the corresponding single and double gluon vertex graphs is shown in Fig. [T5|. As before, 
the linear infrared sensitive pieces in the mass renormalization of the diagrams are identified 
by IR power counting, and will be referred simply as "linear mass pieces" or "linear pieces" 
of the diagrams in the following. Note that the cancellation involves the order g 4 mass shift 
terms from the one-loop diagrams, as will be discussed below. 

Let us first look at the first two graphs of Figs. IT3, and |i"5|b). The vertex of 



|l5|a) is proportional to — ft) (we will always take k the momentum of the first outgoing 
gluon). As discussed for the 1-loop case, the linear mass pieces of the graphs are extracted 
by expanding the vertex around mass shell up to linear terms in the integration variables, 
identified by power counting. In the case of |T5]a) 



CW-W ,. =C5mt(tf-k , -m b )+Cl0mi(tf-k'-m b y (61) 

mass, linear 

The first term in the vertex expansion, Cm\, gives the wave function renormalization piece 
of the diagram. The remaining terms in the vertex expansion contribute to the mass renor- 
malization parts of the diagram. These are at least linear in the offshellness [j> — ft — m b ) , 
canceling a propagator of the graph with an extra factor of i left. This cancellation is as was 
mentioned earlier for the one loop case denoted by a slash on the corresponding propagator 
in the figures. 

The vertex of |l5|b) is a single gluon vertex (Fig. |3|), its expansion up to terms giving 
linear infrared sensitivity in the integration 

\^rtex\ a a \ lmear = Cg m\T« \(j - ftf la + (# - ftf la {$ -ft- ft 1 )+ 

,4 



+ la{P 'ft -01 



linear 



5 1 

C-m A b T a ig la + C-m\T\i -ft- m b )ig la (62) 



Substituting the linearized expressions for the vertices, we find that |l5|a) and pT5|b) cancel. 
The same type of cancellation applies for the next pair of graphs, |T5jc) and |T5]d) 

15a) + [LJ&) = |T3]c) + f[Jd) = 0. (63) 



The case of [15|e-f) and |l5|g-h) is a little more complicated. The leading effective vertex 
is expanded the same way as in fl6l"D , giving 



15e) = Chimt (-zS 22 ) uu - ClOm^ ( C 2 F - -^ > ) uig^ig^ x (64) 



i^m b V 2 



i —i 



d kd kl u u u u Wtt7T m 

f — ft — ft i — m f — ft i — m k A ki 
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Here however the expansion of the single gluon vertex to linear IR sensitive terms in the 
integration results in 3 terms: 



[vertex] 



linear 



CgT b (tf-^^ + ^-k*) 3 ^^- 



7/J - jf - ft)' 
+C—mli 91 pT b 

From (|64l) we get for |l5|f ) 

15/) = -Chim\ (-/T 

d 4 kd 4 k\ 



,10 



(65) 



linear 



C-mti gi pT b + C—mi - - m h ) ig lp T b 



f>=m b 

i 



uu + ClOml Cl - 



CfCa\ _ 



a- -B 



uig^i igy x 



(66) 



-i —i 



—r, — u 1 — u 

— If — If i — m f — f i — m 



ClOmf a 



C F C 



uig^j 



d kd k\- 



x 



— tfi — m^ 13 k 2 kj 



-i —i 



-u 



Unlike the first two pairs of graphs, [l5|e) and [T^f ) do not cancel. The sum of the two graphs 
is the last term of |T^f). Calculating for |T5]g) and |HS|h) in the same way, we find that 



15g) = [THe) 



(67) 



15]e)+|5lf)+|^)+|fi) = 2xC10mUC F 



CfCa\ _ 



uig^i x 



d 4 kd 4 k\- 



-i —i 



m 



,\—m 



l 9iP~gi-gi u 



Next consider the diagrams |T^i)-l). Expanding the vertices around mass shell and keeping 
linear infrared sensitive terms in the integration gives 



15?) = C5im 4 (— zE 2 i) 



t>=m b 



UU 



(68) 



(Note that further expansion of the leading effective vertex would give a quadratically in- 
frared sensitive piece by IR power counting.) 

The single gluon vertex of |IB]j) is expanded to linear accuracy as 



[™rtex\ a a \ Unear = CgT a [f la + f la (rf -#)+••• + la (jj> - Iff 

C-m\T a ig la + C-mlT a ig la {i> - If 1 - m h ) 
% % 

Substituting the vertex expansion (0), |15]j) becomes 
15j) = -C5imi (-zE 21 ; 



linear 



(69) 



+ClQmlC 2 F ui gi a ig^ I d"kd"k x 



p=m b 
a- J 



UU 



(70) 



j4;„ ^4; 



,\ — m 



-Wl3- 



m 
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A similar calculation for |15[k) gives 
15k) = —C5im\ (— zE 2 i 



MM 

■p=m b 



(71) 



— % — z 



— If — m ^ 1 $ — If — Ifi — m ^ tjJ ^ k 2 k 2 
Diagram ITM) is the mass shift term of the tree level amplitude corresponding to E 21 



15/) = C ufu 



Cml 



uu = C5m g 6 2 iAmm 



Adding the 4 graphs and using that S 2 i = a 2 \ + b 2 \g A X we arrive at 
15i) +07) + |I5|fc) +|lf ) = 

-A , -A - % . 1 



(72) 



(73) 



C10m 3 b C F uigj a ig^ / d 4 kd 4 h 



C10m 3 b C F uig-f a I d 4 kd% 



$ — Jf — Jf\ — m 

% a i 



~1 — % 

-igia——u + 



— Ht — m 



— k 1 — )f\ — m 



}f — m k 2 k\ 
— z — % 

wipigia——u 



k 2 k\ 



The calculation for the next 4 graphs, [Hfrri-p) is similar to 1. Note however that in 
the leading effective vertex expansion of |15|m) we must keep the first 2 terms, unlike in the 
case of |T5]i), resulting in an extra term. Apart from this, the difference is only in the gamma 
matrix and color structure. For the sum we find 



15m) +|T^) +[T5|o) + 

-ClOmt ( Cl - ^ 



(74) 



mg'-f x 



d kd k\ 



$ — If — m 



igipigic 



if — lfi—m k 2 k\ 



-i —i 



■U 



-ClOmi Cl - 



CfCa\ _ 



uig^i x 



d 4 kd 4 k r 



-WT 



—i —i 



-ClOmf Ci - 



— If — m 

CfCa 



uig^ a ig~fP x 



—i —i 



<^kd^k x - ig^ a - ig-ip jj-rjU 

f — If — f i — m f — if, i — m k z kf 

Summarizing the leading and single gluon effective vertex diagrams contribute with 
15a) + ■ ■ ■ +^p) = (75) 

ciOml \clui 91 ^ J dVk ^ _ y _Y _ _ \ _ J ria—u 



+ [Cl 



-i —i 



mgr/ igT / a kd k x - Wl*-, OTtjtt" 

f — If — If i — m f — tfi — m k z kf 



+ C F uig~f a / d?kd*k 



WT 



If — m $ 



— % — % 

-i>gipigia——u 



k 2 k\ 
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{ n2 CfCa\ a f ,4, ,4, i i . . —i—i 

+ [ Cf - —) m91 I d M - i - ™ 91 i m l91al91 ^ U 

( n i C F C A \_. r 4 4 i . . i . o-i-i 1 

+ { Cp ~ —) U%91 J d M fc V ~ I ~ m W91a i - ft - m Zg7 FXH 

6. Double vertex graphs Since we keep infrared sensitive terms to linear accuracy, power 
counting shows that the double gluon vertex need not be expanded any further, and so it is 
simplified to 

[double vertex]^ = C10g 2 m 3 b (y a ^T a T b + lpla T b T a ) (76) 

From the above, 

15r) = -C10rnlC F uig^igj a x (77) 

d 4kd4k i j y— u ^la j 77 «57^ jiu 

f — Ip — Ip i — m f — If i — m Ar /cf 

^ — ^ — ft — m a ]/> — ft — m P fcf 

15s) = -ClOm^Clm^ x (78) 

d4 kd 4 ki-} — 5 — rr—72 WlpWla-fijiu 

f — \p — m f — If — f i — m /c z /cf 

-ClOm^ - ^) ui gi " x 

/ d 4 kd 4 k x - j Wl^j 77—72 WlJ91i3-fi-fi 

J f — f — m f — If — tpi — m k z k{ 



u 



15t) = -C10miC 2 F uig-f a x (79) 

/ d^kd^ki- ^ igiaigiPj 77 i 91 P lili u 

J p — If — m f — lfi—m k A k\ 

-ClOml (c f - ^) x 

/ d 4 kd 4 h- — igjpigjaj — 77 i 91 P -r^lJ u 

J f — If — m f — If, 1 — m k A kf 

The sum of the linear infrared sensitive pieces of the double vertex graphs is therefore 

15r) +|I§s) = (80) 



ClOml \ Cpuig^igi* \ d^kd^k^ 7^ — ^ — igja^ 77 ~ i 9^01^TT u 



_^_ft_ m "* '^_ft_ m "« '^2 k\ 
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d^kd^h u u WTI u u WlpWla tt jj 11 

f — if — m f — If — f i — m k z kf 

( 2 C F C A \_. a f ,4 7 ,4, i . i 

+ [ C * ~ —) m91 I d M % - # - m %91 i -t-h- m W91 ^ U 
+ CfrgrT J d^d% J -j-—i giaigi ^ _ ± - m i9^~u 

+ (cl - van* J d*kd% ^ _ * _ - igjpigi y _ £ _ 

Adding ( [75|) and (jSO| ) the sum of the leading, single, and double gluon effective vertex 
diagrams gives 



15a) + • • • + |l5|fc) = (81) 

d 4 kd 4 k ^ - - — igiaigip ^ ZTjj~Z m i9lP ~tf~j$ u = 
-ClOml (-iE^fuu 
Using power counting, the linear term of the sum is (with the use of equation fl58|)) 

15a) + • • • +|T^) = C20rh 3 a 1 b 1 g 4 Xuu (82) 

c. Mass-shift terms from 1-loop graphs The 1-loop diagrams (Fig. |6|) also contribute to 
the order g A cancellations via their mass shift terms. Using the results of section IV A , the 
linear mass renormalization pieces of the 1-loop graphs to order g 4 are 

6a) +§,) +|c) = -C5mt E^uu^ (83) 
= — C 5 (fh + aig 2 + big 2 X^j {aig 2 + big 2 X^ w 
= —Cbm A big 2 Xuu — C40fh 3 aibig 4 Xuu 

( 4 4 \ 5 

6a") = Cm\uu = C fh + a±g 2 + b\g 2 X + a2ig 4 + b2ig A X 

linear \ — ' — ' / 



linear 



UU 



(84) 



linear 



Chm 4L big 2 Xuu + C5m 4 ^i*? Aim + C20m 3 aibig Xuu 



,2 



As we have shown before, the order g terms of the 1-loop diagrams cancel. The mass shift 
terms due to the order g 4 self-energy corrections were already included in the calculation. 
The order g A terms we have not included so far are the "cross terms" , those proportional to 
axb\. These are the terms we denoted with the expression "mass shift terms" on Fig. fT5jx), 
summing up to 

15x) = -C2$nV i a 1 b 1 g A Xuu (85) 

This contribution cancels the remaining linear term of the 2-loop vertex diagrams of group- 1 
and group-2, (|8"2"D, so the linear infrared sensitive terms in the mass renormalization of the 
first two group of diagrams with the mass shift terms from the 1-loop diagrams cancel 

15o) + • • • = 0. (86) 
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2. Group-3 and the group of 1-particle reducible diagrams 



The linear mass renormalization parts of the group-3 diagrams cancel with the corre- 



sponding pieces of 1-particle reducible graphs, as shown in Fig. 16 



Expanding the leading effective vertex of |l6|a) as in (|6l|) , and the single gluon vertex of 
|16|b) as in (|69|) to linear accuracy it is easy to see that [T6|a) and |T6|b) cancel. The next two 
pair of graphs [l6|c) and |l6|d) cancel similarly: 



16a) + = |ljc) +[T6]d) = (87) 

The mass shift term of the tree amplitude, corresponding to £ 2 3, |l6|g) 

16g) = Cm\ uu = ^>m% 2 z9 ^uu (88) 

cancels with the remainder of the one-particle reducible group, as we will show below. For the 
1-particle reducible diagrams the leading effective vertex of |l6|e) , |l6|f) and the corresponding 
mass counterterm graphs is expanded to first order in the offshellness 

(i-tyf =C 5mt^-^-m b ) (89) 

mass, linear 

Keeping the second term as in (|61|) would result in a quadratic infrared sensitive term by 
IR power counting. From the vertex expansion (|5^) we get for the linear mass pieces of the 
diagrams 



16e) +H/0 + • • • = -Cbimi (Z- 1 - l)^ (-iEi)^ uu (90) 

The ellipsis in the formula denote the inclusion of the mass counterterm diagrams, as in Fig. 
16. Here we used the one- loop Ward-like identity (Fig. |5[) and the perturbative expansion 
afZ^.Eq. ©. 

For the remaining one-particle reducible diagrams |16|h-k) the expansion of the single 
gluon vertex up to linear infrared sensitive terms in the integration gives one term 

[vertex] a a \ Unear = C 5 -m^ Q T a (91) 

t 

Substituting the vertex expansion (0), we find that 

I6/1) + 00 + • • • = Hlb') + • • • = C ^ m t i 2 ' 1 - l ) g2 i-^i)\p= mb uu. (92) 

Thus the sum of the linear pieces of the 1-particle reducible graphs is 

16e) + • • -|l|fc) + • • • = Chim\ (Z' 1 - l) 2 (-zSi 
= C5m 4 (Z- 1 - l) b ig 2 \uu 



uu (93) 

ft=mi,,linear 



The sum (|93|) cancels the leftover S 2 3 mass shift piece, (^8|), since by the one loop Ward-like 
identity (Fig. R[) 
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-h (z- 1 - 1 



and Eq. (|8BD (Fig. $Bg)) becomes 

16g) = C5m 4 6 23 (y' 4 AMn = —Chfn%ig 2 (z~ x — 2 Xuu, 



(94) 



(95) 



Hence all linear mass renormalization pieces of the group-3 and the 1-particle reducible 
diagrams cancel: 



16#)+ye) + ---|Jfc) + 







(96) 



3. Group-4 

The cancellation within group-4 (Fig. [T7]) proceeds in the same way as for the 1-loop 
diagrams (see Fig. || and section [IV A| ). 



4- Group- 5 



Fig. [18| shows the cancellation between the group-5 diagrams. We can calculate the linear 
parts of the first 4 diagrams |l8|a-d) expanding the leading effective vertex to second order 
and the single gluon vertex to first order in the offshellness. A straightforward calculation 



gives 



18a) + 



0c) +[TB]d) = 



(97) 



The case of the next 4 diagrams, |18|e-h) is more interesting. |18e), the mass-shift piece of 
the tree amplitude corresponding to £25 is 



18e) 



Cm\ 



uu 



5Cm &255 1 Xuu 



(98) 



Denoting the triple gluon vertex with S^ s , the next diagram, fLqf) gives 
18/) = uig-i a T a J d A k b d A k c - 



(99) 



[vertex] 



5rj-\c nobc 



-iir\°T c 



, ,,, m- ,6 k 2 b h\ (/.v. • k c f 

Sa% = [9*13 (ka - h) S + 9I3S (h - k c ) a + g Sa (k c ~ K) p 

The single gluon vertex is expanded to linear accuracy as follows: 



[vertex] 



3.b 
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linear 



C-m A b i gi p T b + C— 

+C—mlig^T b 

i 



in,. 



-fa- If 



, m b ) ig^T b 



(100) 
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Substituting the linearized vertex, |18|f) gives 3 terms 
18/) = -CUmi HS 25 )|^ W6 uu + ClQm 3 ui gi a T a ig^T b x 



;ioi) 



d A kbd A k c - 



- H-c. — m 



-*gYT c s, 



-i —i 



^ (fib ^c) 



Cl0m 3 uig-? a T a / d 4 k b d 4 k, 



- #b - tyc - m 



ig^T h i gi 5 T c S. 



abc 



—I —I 



^ ^b (^b ~t~ k c 



To get the linear terms of I8"|g) and p!8[h), we need to keep only the first term in the 
expansion of the double gluon vertex. Thus 



18g) = C10m 3 b g 2 u \ 1 a 1 f3 T a T b + 1 P 1 <*T b T a 



d 4 k b d 4 k c 



Srpc oabc 



m 



-igi 6 T c s, 



x 

i —i 



(102) 



al3S e b ki (k b + k c y 



u 



and 



18h) = C10m 3 g 2 uig~f a T a x 



d 4 k b d 4 k c - 



(103) 



s: 



abc 



— I —I —I 



- # 6 - # c - m 

Using that £ 2 5 = a 2 ^g 4 + b 2 ^g 4 \ ) the sum of the group-5 graphs gives 
18a) + ---+|T|/i) = 



(104) 



-C10m^^ 7 a T a / d 4 k b d 4 k 



-i —i 



-Cl0m 3 uig^T b ig-? a T a / d 4 k b d 4 k, 



-i — i — z 



&rpc nabc 



- ktr.-m 



-i gi °T c S t 



^ k"b fob kc) 



This remainder can be shown to be zero (up to linear accuracy) by renaming variables for 
example in the first integral for k' c = (k b + k c ) and k' b = —k c . Therefore the linear mass 
renormalization pieces of group-5 diagrams cancel. 



V. SUMMARY AND CONCLUSIONS. 

In this paper we have shown by explicit perturbative calculations upto the second loop 
order, that when the inclusive decay width of a heavy quark is expressed in terms of the short 
distance mass, there is complete accord with the operator product expansion and the leading 
power corrections are of order 1/m 2 . This is the first such calculation to the second loop 
order that checks the dependence on the infrared cutoff to linear accuracy for the non-abelian 
theory. Because of the non-abelian interactions there could have been potential problems 
due to the singularity in the t channel for processes involving the forward scattering of soft 
gluons. However we see explicitly that these do not hinder the cancellation of the infrared 
sensitive terms to linear accuracy. 

As emphasized earlier, the above result is true not just for the semileptonic decay con- 
sidered in this paper but for any inclusive decay, like for example the radiative one. Such 
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a cancellation was linked to the KLN theorem which had been suggested earlier to be 
the general principle behind this for inclusive enough observables. Since for the inclusive 
decay, there is only a single particle in the initial state, our calculations enjoyed several sim- 
plifications. It is known (Kj that for at least two colored particles in the initial state there 



is a breakdown of the Bloch Nordsieck mechanism in perturbation theory. Thus it would 
be interesting to investigate the fate of the 1/Q corrections for the Drell Yan process. For 
this process the initial state is no longer trivial and in addition furthur new complications 
arise in higher loops due to the mixing of soft and collinear singularities ||11| . The 1/Q 
corrections are known [T^j to cancel to one loop order and the situation at higher loops is 



an open question. The techniques developed in this paper would be useful for this study as 
well, and we hope to return to this question in a future publication. 
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APPENDIX A: CANCELLATION OF INFRARED SENSITIVE TERMS TO 
LINEAR ACCURACY IN THE WAVE FUNCTION RENORMALIZATION: 

GROUP-5 



For the last group of diagrams (Fig. [12]) the presence of the non-abelian vertex makes 
the algebra a little more complicated. The calculation proceeds essentially along the same 
lines as for the second group of diagrams (see section |111 Bp . We would like to show that 
for the logarithmically and linearly infrared sensitive parts (referred to as "logarithmic and 
linear pieces" in the following) of the diagrams of Fig. [T2| 



[12a) +W>)La V efuncUon = ^1 [Z' 1 - l) ^ UU (Al) 

= — |l2]c) + ^fi)] wave f unction 



As before, the logarithmically and linearly infrared sensitive pieces of the integrals are iso- 
lated by IR power counting. In the rest frame of the heavy quark p = (m b , 0, 0, 0), and 



u^ 



u u^ 



u 



f z -l _ ] \ = °P» t>=m b = 9p ^ mb 

V /9 4 ,S25 vrfu uu 

so we would like to show that for the logarithmic and linear terms in the integration 



i>=m b 



dp, 



o 



i>=m b 



U = -[12c) +mi)Lavef unction ( A 3) 



As before, the derivation of £25 with respect to po results in an insertion of an extra vertex 
27° and an extra propagator in all possible way. The resulting diagrams are similar to 
|T2"|c) and |T2|d) with the Cm\ vertex replaced by — Cm^7° (shown in Fig. [TB| ). Writing out 
explicitly 
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Cmlu 



as 



25 



dp 



u 



5 C F C A 
-Cm b — • 



j>=m b 
% 



(igf x / d^fad^fa^x 



(A4) 



-2p (fa + fa) + (fa + k c ] 
i 

-2p (fa + fa) + (fa + fa)' 



-i —i 



-2pk c + kl kl k\ (fa + faf 

2 

—I —I —I 



+ / d A fad 4 k c N 2 x 



-2pk c + k 2 c 



k 2 k^ [kb + k c ) z 



Denoting the triple gluon vertex by S a ps (where the color factor and the coupling are taken 
out) and k a = — (fa + k c ) 



Nt = (p' + ^ a + m) 7° (p 1 + # n + m) ^ (J - # c + m) ^ s S aPS u 
N 2 = u-f a (p l + l^ a + m) ^ ($-ty c + m) 7° (p 1 - ^ c + m) -f S S af35 u 
Saf3S = 9a/3 (k a - fa) s + gp S (fa - fa) a + g Sa (fa - fa)p 

Comparing this to diagrams [L^c) and |12]d) 



(A5) 



-12c) = -Cm\ 



CfCa 



(igf / d%d 4 faN lc x 



(A6) 



-2p(fa + fa) + (fa + fa) _ 



-i —i 



2pfa + kl kl k 2 b (fa + fa) 2 
(igf j d 4 fad A k c N 2d x 

2 



-2p (fa + fa) + (fa + faf [-2pfa + k 2 c 
with the numerator factors 



-i —i 



kl k% (fa + k c 



N lc = u^ a (p' + ^ a + m) (p' + ^ a + m) ^ (p'-^ c + m) i & S aps u 
N 2d = uj a (p' + ^ a + m) ^ (p'-k' c + m) (p'-^ c + m) j s S a p S u 



(A7) 



we have to show Aq = Aq c and N 2 = N 2d to linear order in the integration. By power 
counting we must keep terms to the order k in the numerators in addition to the triple gluon 
vertex k dependence. The logarithmic terms arise dropping the k dependence everywhere 
in the numerators except in the triple gluon vertex S a p$. A simple calculation shows that 
the numerator factors for the logarithmic terms are zero. For the linear terms, a little more 
calculation gives 



N\ = Am 2 {(fa - fa) Q utfaU + 2 (fa - fa) ulf, c u+ 

+U (ty b - # c ) IjtaU + Utya (# c - # ) U - uk\ c (tf a ~ ty b ) U - U 



(A8) 



l a ) ki c u} = N : 



lc 



where the rest-frame condition was used. An entirely similar calculation gives N 2 = N 2( i to 
linear accuracy. Thus the numerator factors are equal up to linear order in the integration, 
and equality ( |A3[ ) is established to linear order. Hence we have shown the cancellation of 
the logarithmic and linear wave function renormalization pieces of group-5. 
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FIGURES 




FIG. 1. The forward scattering amlitude for the heavy quark decay whose imaginary part 
gives the decay width. Shaded blob denotes the hard part, while S denotes the soft gluon interac- 
tions. 
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FIG. 2. a) Emission of a single gluon from the final state quark. The blob denotes all 
possible gluon self-interactions before coupling to the heavy quark, b) An example for the gluon 
self- interactions included in a), c) Two gluon emission from the final state quark. 
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a, a /3,b 
c) 

FIG. 3. a) leading effective vertex, b) single gluon effective vertex , c) double gluon effective 
vertex 
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a) b) c) 

FIG. 4. Wave function renormalization and Bremsstrahlung diagrams to order g 2 . 




FIG. 5. One-loop order Ward-like identity 




FIG. 6. Mass renormalization and the corresponding bremsstrahlung diagram to order g 2 . 
Figure a) depicts the mass renormalization piece of the graph. The propagator with slash is 
cancelled by the expansion of the leading effective vertex. The subscript Si of d) denotes the mass 
shift due to the order g 2 self-energy correction. 
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FIG. 7. Self-energy corrections of order g 2 and g 4 
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FIG. 8. The first group of 2-loop diagrams corresponding to £21. 
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FIG. 9. The second group of 2-loop diagrams corresponding to £22- 
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d) 



FIG. 10. The third group of 2-loop diagrams corresponding to £ 
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b) 



FIG. 11. The 4th group of 2-loop diagrams corresponding to £24- The blob denotes fermion, 
duon and ghost loops. 
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b) 



d) 



FIG. 12. The 5th group of 2-loop diagrams corresponding to £ 



25- 
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FIG. 14. 1-particle reducible diagrams for wave function renormalization. The ellipsis denote 
all possible inclusion of mass counterterms. 
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mass shift terms linear 
/ W= / \ + from 1-loop 

graphs 

x) 

FIG. 15. Mass renormalization of order g 4 : group- 1 and group-2. The linear infrared sensitive 
mass renormalization pieces of the leading effective vertex graphs cancel with the corresponding 
single and double gluon vertex graphs and with g 4 order mass shift terms from the 1-loop diagrams. 
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a) b) c) d) 




FIG. 17. Mass renormalization of order g 4 : cancellation of linear pieces in group-4 
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FIG. 18. Mass renormalization of order g 4 : cancellation of linear terms in group-5. 
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